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Abstract 

Belle data on 77 — > tj'tttt are refitted using a broad J PC = ^ peaking in the mass 
range 2250-2300 and X(1835), but without 7/(1760). There is the possibility that the broad 
h signal may be identified with the ^ glueball predicted originally by Morningstar 
and Peardon. The X(1835) is confirmed to have a resonant phase variation. 

PACS numbers: 14.40. Be, 14.40. Df, 12.39.Mk 

1 Introduction 

The Belle collaboration presents new data on 77 — > ?/7r + 7r~ In their Fig. 3, there is 
evidence for fine structure in the 7/V7r mass range 1700-1900 MeV. On a larger scale, there is 
a conspicuous broad peak centred at 2300 MeV with a full width at half-maximum of ~ 750 
MeV, see their Fig. 2(b). As a shorthand, this peak will be called V(2300). No fine structure 
is visible in this broad peak from 2000 to 2800 MeV. It deserves attention, since it could be the 
0~ glueball predicted by Morningstar and Peardon near this mass [2]. 

It would not be surprising if this glueball is very wide. Zou, Dong and I have drawn attention 
to a very broad J PC = 0~ signal observed in J/ip radiative decays j3]. It accounts in a simple 
way for successive peaks in J/ip — > 7X, where X — > pp, ww, K*K* and 00 channels, with 
flavour-blind coupling strengths. J/ip radiative decays are dominated by 7GG, where G are 
gluons. The conclusion of Ref. [3] was that there is a broad h signal consistent with a 
glueball with mass M = 2190 ± 50 MeV and width T = 650 ± 100 MeV. The half-width of the 
lower side of the peak in Belle data is 300-350 MeV. It would be a pity to miss the 0~ glueball 
if it is really there. 

Belle base their analysis on the claim by DM2 to observe 7/(1760) in data on J/ip — > 
7(7r + 7r~7r + 7r~) [4j. However, an analysis of Mark III data on the same channel showed that 
the 1760 MeV peak has J = ++ , though it does sit on a large, broad 0~ background [5]. 
[A technical detail is that there is no interference between ++ and _+ after summing over 
spin orientations of the J/ip.} Peaks at 1500 and 2105 MeV were also fitted with J p = + . 
Furthermore, high statistics data of the Fermilab E760 experiment [6] on pp — > (7777)71"° fit all 
three peaks accurately with the same mass and width for these J p = + states; J p = 0~ is 
forbidden in rjr) by the Pauli principle. Many authors have been confused by the fact that the 
PDG [7] does not mention Ref. jl] under 7/(1760), though it is listed under / (2100). 

The existence of 7/(1760) therefore rests on (a) the BES I analysis |8J, where M = 1760 ± 35 
MeV, T ~ 250 MeV, (but quoted as not well determined), (b) the BES II analysis of J/ijj — > 
7(wu), M = 1744 ± W(stat) ± I5(syst) MeV, T = 2U±H ± 25 MeV |9 . There is a serious 
objection to this second source. A well known relation, coming from SU(2) symmetry, is that an 
isospin 1 = resonance should have equal couplings to uu and p°p°, because light quarks do not 
discriminate between charges. There are three charge states for pp, so the relation is normally 



written g 2 (pp) = 3g 2 (u)u), where g are coupling constants. This relation applies equally well to 
qq states, hybrids and glueballs, which all obey SU(2). 

The branching fraction quoted for production of ?y(1760) in the BES II 'yuu data is (1.98 ± 
0.08(stat) ± 0.32(syst)) x 10~ 3 . This is larger than the p°p° weighted mean branching fraction 
from DM2 and Mark III over the entire mass range up to 2 GeV, namely ~ (1.23 ±0.25) x 10~ 3 . 
It should lead to a huge pp peak at 1760 MeV, in disagreement with BES I, DM2 and Mark 
III data. The DM2 collaboration did claim a small ?y(1760) signal in data on pp — > 7r + 7r~7r + 7r~, 
but without observing any phase variation. Their branching fraction was a factor 4 smaller 
than BES II claim in 'yuu. The 7r + 7r~7r + 7r~ data are experimentally much cleaner than 'yuu, 
u — > 7T + 7r~7r°, where there are 5 photons, hence large combinatoric problems. 
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Figure 1: Data and fits to J/ip — >■ (a) tjtctc, (b) pp, (c) wu, (d) K*K*, (e) (jxp and (f) all 
channels, from Ref . {3j . Points with errors show averages of data from BES I, DM2 and Mark 
III. Reproduced with permission. 



Fig. 1 reproduces the magnitudes of branching fractions used in Ref. [3] for five channels and 
the total. This analysis used fully analytic amplitudes where the denominator of the amplitude 
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takes the form 



D(s) = M 2 -s-m{s)-iY J 9 2 3 P 3 ^) (!) 
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S -M 2 r°° Ejgjpji/) ds' 

m{s) = —^ p Lr (P - s)(s> - m (2) 

Here P stands for the Principal Value integral; sthr is the threshold for each channel j. Note 
that this is not an 'optional extra'; it is a requirement of analyticity. Wherever an opening 
channel produces a peak, there is a cusp, which need not be resonant. 

The pp channel shown in panel (b) peaks just below 1600 MeV because of the L = 1 centrifugal 
barrier for production. In panels (b) and (c), amplitudes are restricted to obey the SU(2) 
relation. The peak in ww can be confused with a resonance, but at half-height it has a width of 
~ 350 MeV. In Ref. [3] it was shown that there was no pole in the ujuj amplitude, nor in i]tttt in 
this mass range. 

Section 2 discusses first a fit to the broad peak considering only J p = 0~ — > rja. The objective 
is to refit the Belle data without ^(1760). The X(1835) component is needed to fit the mass 
distribution. Note that the [770"] decay (where L is orbital angular momentum) is forbidden 
for 77—7- 1 ++ (Yang's theorem). 

However, this is not the whole story. Belle's Fig. 6 displays the 7t + tt~ mass distribution 
for rj'im events in the mass range 2200-2700 MeV/c 2 . There is an / 2 (1270) peak in 77?t?t, 
with an intensity ~ 85% of that of the a near 1 GeV. A likely amplitude producing the fi is 
y(2300) — > [r]' f^\L=2, with the result that there is a single broad 0~ resonance. This possibility 
is investigated in sub-section 2.1. 



2 Initial fit to the r) f ir + 7r peak 

Ideally, one would adopt the same approach as used in Ref. [3], where a coupled channel fit to 
many open channels was made. However, presently there is limited knowledge of some important 
amplitudes in the mass range above 1800 MeV, e.g. 77 — > KKtt and t/tttt. Only the simplest 
parametrisation for the broad component visible in Belle data can be used at the moment, but it 
is essential to accomodate the opening of rj'a and r]' f 2 phase space p. The simplest Breit-Wigner 
amplitude f(s) with these features is 

f(s) = MT{s)/[M 2 -s-iMT ] (3) 
r(s) = FF(k)g 2 (7]'mr)p(7] , Trn). (4) 

If other decay modes of V(2300) exist, such as KKn, 777171", uw and pp, then To will be summed 
over all decay modes. Here it is taken as a constant, the simplest possibility. For the numerator, 
a form factor FF(k) = exp(—ak 2 ) is used, where k is the momentum of the 77' in the overall 
centre of mass. Good fits are obtained for a = 1.5 to 3.0 (GeV/c) -2 , consistent with data on J/ip 
radiative decays [3]. The form factor arises from convolution of form factors for the outgoing 
777T7T final state and the initial 77 interaction. The phase space factor p is obtained from the 
integral over 3-body phase space, given by equations (39.19) and (39.20) of the Particle Data 
Book [7]. For fits where the 0~ initial state decays to both r/'a and vj ' f%, followed by decays of 
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both channels to 7r + 7r~, T(s) of Eq. (3) needs to include fully coherent interferences between 
both channels. The a is parametrised by the amplitude given in Ref. [TO], Table 1, entry (iii); 
this parametrisation allows not only for decays a = /o(500) — > nn, but also to KK (which is 
quite significant), 7777 (small) and Arc (large above 1350 MeV, but affecting only r)'7nr masses 
above ~ 2250 MeV and rather uncertain in magnitude). The X(1835) is included in the fit 
multiplied by an isobar model phase factor exp(i5), with constant 5. 




mass (GeV) 



Figure 2: Fits to Belle data for 77 — > ri'7i + -n~ from their Fig. 2(b); full line: fit with F(2300) 
and X(1835); dashed curve, the fitted intensity of X(1835); dotted curve, fit with only V(2300); 
chain curve, the full intensity including [r]' f?\i,=2- 



Three curves on Fig. 2 show the fit to Belle data, initially using only 0~ — V rja over the 
whole mass range up to 2800 MeV; a fourth curve shows the effect of including 0~ — > [r//2]L=2 
as described in sub-section 2.2. The full line shows the fit including V(1835). Its mass and 
width are allowed to vary within the statistical and systematic errors quoted by BES III 
A resonant phase variation is required for V(1835) to account for interference with V(2300). 
The mean x 2 P er point is 1.25 after allowing for the 7 fitted parameters. Two points at 1.69 
and 1.73 lie 2.8 and 2.4 standard deviations above the fit, but would require a very narrow peak 
inconsistent with the broad ww cusp. The dotted curve shows the fit without V(1835). The 
contribution from V(1835) to the full line is an 8.2 standard deviation effect, after correcting 
deviations from the fit by dividing by 1.25 and allowing for the change in the number of fitting 
parameters from 7 to 4. So the X(1835) is confirmed as a 0~ + resonance. 

Using the whole mass range, the mass fitted to the broad 0~ signal in Fig. 2(a) is 2300^|o 
MeV and T = 750^ MeV, where the signs display the correlation between M and T. The 
dashed curve at the bottom of Fig. 2 shows the optimum contribution from V(1835). The 
r/(2300) intensity without V(1835) is shown by the dotted curve. 

If 77 couples to wu, as seems likely, the full width of the cusp due to this threshold is ~ 350 
MeV at half-maximum from Fig. 1(c). This will alter the entire 7/(2300) mass distribution 
slightly, but cannot be predicted without the 77 — > ww coupling constant. 
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2.1 The effect of rf f 2 (!270) 

The J p = 0~ — > [rj'f 2 ]L=2 contribution provides a ready explanation of f 2 production. Their 
coherent sum is deduced using the relative magnitudes of f 2 and o signals in Fig. 5 of Belle. 
The phase space factor for f 2 includes a standard L = 2 centrifugal barrier factor B(s) = 
9/c 4 /(9 + 3k 2 R 2 + k A R A ). The value R = 0.725 fm is used. This assumes the same radius of 
interaction as the Gaussian form factor FF = exp—[k(GeV/c)/hc] 2 R 2 /Q = exp(— a/c 2 ), using 
the optimum a = 2.25 (GeV/c)~ 2 . Following the isobar model approach, the amplitude requires 
a factor exp i5, where 5 is a constant, 25°. The angular distribution between rf and f 2 is isotropic, 
like rfa. A summation is made over the f 2 decay. 

The \rf f^L=2 amplitude produces a slow rise centred at 2400 MeV. In order to accomodate 
this rise, the width fitted to F(2300) decreases by ~ 100 MeV to 650 MeV, with an associated 
reduction in the mass of ~ 50 MeV. These values provide an estimate of systematic errors for 
the mass and width of V(2300). This extra component allows slightly more freedom in the low 
mass range considered in the previous section. The significance of the X(1835) contribution falls 
from 8.2 to 7.1 standard deviations, but there is no essential change to the fit, which is shown 
by the chain curve in Fig. 2. It produces (i) a slightly larger X(1835) peak, (ii) a V(2300) mass 
reduced to ~ 2250 MeV, and (iii) a higher tail near 2700 MeV caused by the rising centrifugal 
barrier. It does produce an angular distribution resembling Belle's Fig. 5. However, a full fit 
to the Dalitz plot will be required to be precise about the angle and energy dependence of this 
term, if indeed it is present. 

Belle suggest a contribution to the V(2300) peak with J p = 2 + . This would produce the final 
state T7'/ 2 (1270) L=1 . There is also the possibility of J p = 2" — > [rj'f 2 ]L=o- Only an analysis of the 
Dalitz plot can identify such contributions. A remark is that the angular momentum analysis of 
these cases needs to obey gauge invariance for the photons. In their centre of mass, they have 
only helicity amplitudes |1, 1 > and |1, — 1 >, but no |1, > component. If axes are used in the 
77 rest frame with the z-axis along the direction of the photons, only the x and y components 
of the spin 2 combination contribute; their intensities add incoherently. To describe the f 2 it is 
necessary to rotate axes to its direction and then make a Lorentz boost to its rest frame. Using 
rotation matrices and Clebsch-Gordan coefficients is prone to mistakes, but possible. A simpler 
approach is the so-called method of Wick rotations. Axes are first rotated to the direction of the 
f 2 using angles 9, (f>. After expressing momenta of pions from the f 2 decay in this system, the 
Lorentz transformation of the pions is made to the f 2 rest frame. Finally, the axes in that frame 
are rotated back through angles —9, —<p, taking care that the product of the two rotations is the 
3x3 matrix with unit diagonal elements. The two rotation matrices cancel. The Lorentz boost 
changes the angles of the pions from the f 2 decay between the 77 rest frame and the f 2 rest 
frame, but the usual Clebsch-Gordan coefficients and spherical harmonics give the correct decay 
amplitude, which also requires the usual centrifugal barrier factors. A check on the procedure is 
that all amplitudes, including those for rf-n-n used above, are orthogonal when integrated over all 
angles. If the experimental acceptance is included using the Monte Carlo simulation, the effects 
of cuts and acceptance are immediately apparent in the interferences between amplitudes. The 
third method is to use covariant tensor expressions for amplitudes, but this is unfamiliar to most 
experimentalists. 
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3 Concluding Remarks 



The glueball hypothesis is clearly a matter of conjecture, but is worth following up with further 
studies of 77 — > rjnn and KKtt in particular. These are important for checking whether the 
V(2300) decays flavour blind. The observed broad, well defined peak in ti'-k-k does not look like 
a conventional qq state, where masses are typically below 300 MeV. The lower side of this peak 
requires interference with X(1835) with J PC = 0~ + in agreement with BES III; the X(1835) 
has a definite resonant phase variation. 

At present, there is no clear benchmark for the mass scale of the 0~ glueball, but Lattice 
estimates are in the mass range 2250 [12] to 2590 MeV [2j. Mass ratios are better predicted than 
absolute values. Morningstar and Peardon predict that it will have a mass 1.50 ± 0.04 times 
that of the + glueball and 1.08 ± 0.04 times that of the 2 ++ glueball. A second + glueball is 
also predicted in this mass range. It is possible that the 2 + glueball is to be identified with the 
/2(1950), which is observed in decays to 7777, 4-7T, KK and KKirn; it has a width of 472± 18 MeV 
[7j. There is a new quenched Lattice Gauge calculation of the glueball spectrum by Gregory 
et al. [13], who give a careful review. The present theoretical situation is that eigenvalues, 
i.e. masses, have been calculated. These are based on couplings to gluons. What is presently 
not clear is the effect of decays to qq nonets. These are in most cases still buried in two point 
correlation functions which include glueball and qq combinations [T4"] . 

In BES III data for the 7/ channel, there are further peaks at 2122 and 2376 MeV with 
widths of 83 MeV for both. A natural interpretation of the A(1835) is that it is the n = 3 
radial excitation of 7/(958) (n = 1) and 77(1405/1475) (n = 2). For the latter, the evidence for 
two separate states is not conclusive; the stronger 77(1475) decays dominantly to i^i^* (890) in a 
P-wave, and the k 3 increase of the P-wave shifts the average mass of the peak up by ~ 35 MeV. 
Also Wu et al. [15] have proposed an interpretation of 7/(1405) in terms of a triangle graph where 
KK from KK* decay rescatter to rjTf. In J/tp — > 7(77^7?), there is a large dispersive peak in Fig. 
1(a) at 1500 MeV producing a strong enhancement of 7/(1405) — > tj^tt. Achasov and Shestakov 
suggested in 1985 a natural explanation for a very broad J PC = 0" signal observed in J/ijj 
radiative decays to pp and uu [17]. In a later paper they suggested why 77(1440) and 7/(1475) 
are not observed in 77 collisions [15] ; further consequences of this suggestion were studied there. 
A recent third paper reviews the question comprehensively including the latest data and makes 
recommendations for further work [T9~] . 

It is not yet established that the peaks at 2122 and 2376 MeV have J p = 0~. However, if that 
is the case, mixing with the broad Y (2300) would enhance their visibility. The gluon interaction 
is likely to be of short range, judging by the funnel potential, but mixing with qq components 
which peak at larger radii reduces the zero-point energy. The sequence of peaks from r/(958) 
to 2376 MeV lies close to a straight trajectory of M 2 v 77, with a slope of 1.18 GeV 2 , like that 
observed in Crystal Barrel data for many resonances, namely 1.143 ± 0.013 GeV 2 [16J. 

In summary, the Belle data can be fitted well with just X(1835) and a broad J PC = h signal. 
An f2 component arises naturally from 0~ — > [f]'f2]L=2 but that hypothesis needs confirmation 
from a full analysis of the Dalitz plot. 
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